• 专利附录
  • 专利说明
  • 权利要求
  • 类似技术
  • 同族专利
  • 引用文献
  • 法律状态
  • 优先权
    • 专利摘要:
    A method for generating a seismic image, comprising the steps of providing an input trace corresponding to seismic signals received at a receiver at a known receiver depth, using the input trace to generate a series of angle traces, each having a known time dip, for each angle trace, determining a filter
    公开号:AU2013215333A1
    申请号:U2013215333
    申请日:2013-01-29
    公开日:2014-07-31
    发明作者:Jonathan Lewis Sheiman；Ailin YIN
    申请人:Shell Internationale Research Maatschappij BV；
    IPC主号:G01V1-38
    专利说明:
    WO 2013/116195 PCT/US2013/023583 ANGLE-DEPENDENT DEGHOSTING OF SEISMIC DATA RELATED CASES Not applicable. FIELD OF THE INVENTION [0001] The invention relates to a method for deghosting seismic data acquired by a marine seismic source and receiver assembly, and more particularly to a method in which the effects 5 of seismic reflections by the water surface, or ghosts, are mitigated with a angle-dependent deghosting algorithm. BACKGROUND OF THE INVENTION [0002] Marine seismic data is typically acquired by hydrophones in a cable or cable array that is 10 towed behind the acquisition vessel and kept several meters below the sea surface. The hydrophones record both waves which travel directly from the earth to the receivers (the "primary" signal) and those same waves after they have been reflected at the sea surface and arrive after a time delay (the "ghost" signal). The ghost signal interferes with the primary signal constructively or destructively, giving a distorted signal. This is particularly true when the source and receiver are relatively close together, as 15 the time delay of the ghost is not sufficient to separate it from the primary signal. Ghosts can occur relative to the source ("source ghosts") and relative to the receiver(s) ("receiver ghosts"). [0003] Because ghosts distort the primary signal, they reduce the value of the resulting seismic data. Therefore, it is common to undertake various techniques to reduce the effect of ghosts. Mitigation or elimination of ghosts from the desired data is called "deghosting." 20 [0004] If data are gathered simultaneously using both hydrophones (pressure sensors) and geophones (accelerometers), the difference between the two datasets can be used to deghost the data. If two types of sensors are not available, however, other deghosting techniques must be applied. [0005] Some conventional seismic processing techniques attempt to deghost the data using an assumption that the waves are traveling in a near vertical direction as they approach the instruments. 25 However, waves can and do travel at angles through the water, due to reflector dip and/or shot-receiver offsets. The nature of the distortion varies, and the vertical deghosting procedure doesn't fully remove it. The problem becomes more significant as the depth of the receivers increases. Thus, it is increasingly desirable to provide a deghosting technique that takes into account the propagation angle of the incoming seismic signal. 30 WO 2013/116195 PCT/US2013/023583 SUMMARY OF THE INVENTION [0006] In accordance with preferred embodiments of the invention there is provided a deghosting technique that takes into account the propagation angle of the incoming seismic signal. [0007] In one preferred embodiment, at least one input trace corresponding to seismic signals 5 received at a receiver at a known receiver depth is provided. The input trace is used to generate a series of angle traces, each having a known time dip (px, py) = (dT/dx, dT/dy). F1 2 2- 2) Then, for each angle trace, a filter H is determined using sin(O) = 1 - v p +p where v, is the propagation velocity in water to determine an angle 0 and H (co) = 2i sin(iaoz, cos(6) /vv)where z, is the receiver depth and ow is the angular frequency. 10 Each angle trace is then corrected for ghosting by using the inverse of the filter H so as to generate a plurality of deghosted angle traces. Each deghosted angle trace is then copied over to a plurality of output locations, applying in each case the time dip for that angle trace, so as to generate a plurality of corrected angle traces and, for each output location, a plurality of the corrected angle traces at that output location are summed so as to generate a replacement 15 trace. The replacement trace can be used to generate a seismic image. Alternatively, steps a) - g) may be repeated for a plurality of input traces, in which case step g) may comprise, at each output location, summing the replacement traces at that output location so as to generate an output trace and using the generated output traces to generate a seismic image. [0008] The output locations may or may not correspond to the input locations. The number of 20 input traces may or may not be greater than 100. Step b) may or may not include generating at least 30 angle traces from each input trace. [0009] Steps a) - i) may be repeated after a period of time so as to generate a second seismic image and time-spaced seismic images may be compared so as to generate a time-lapse image. [0010] The seismic signals are preferably received at a receiver selected from the group consisting of hydrophones, geophones, fiber optic sensors, ocean bottom sensors, streamers, and combinations thereof and are preferably received at a towed marine receiver. [0011] As used in this specification and claims the following terms shall have the following 25 meanings: "offset" refers to the horizontal distance from the source-position (or from the center of a source group) to a seismic detector, or more commonly to the center of a detector group measured in the horizontal (x-y) plane. For a vertically staggered source array, the source- WO 2013/116195 PCT/US2013/023583 position is assumed to be located at a point in the x-y plane directly above the center-point of the stacked elements; "in-line" refers to a line whose axis lies in the horizontal (x-y) plane and lies parallel to the direction in which the 3D seismic data was acquired; 5 "primary signal" refers to that portion of the sound energy that propagates from the reflector(s) to the receiver(s) without being reflected off the sea surface; and "ghost" refers to a delayed negative representation of the primary signal that is the reflection of the sound energy by the sea surface. 10 BRIEF DESCRIPTION OF THE DRAWINGS [0012] For a more detailed understanding of the invention, reference is made to the accompanying Figures wherein: Figure 1 is a schematic illustration of an exemplary system that could be used to collect marine seismic data for use with the present invention; 15 Figure 2 is a geometric illustration showing the dependency of ghost signals on the angle of propagation; Figure 3 is a schematic illustration of the concepts of the present invention; and Figure 4 is a flow chart schematically illustrating one embodiment of the present methods. 20 DETAILED DESCRIPTION OF A PREFERRED EMBODIMENT [0013] Referring initially to Figure 1, a seismic exploration vessel 10 is shown towing a marine streamer cable 12 through a body of water located above the substrata 20 that are to be seismically explored. Cable 12 can be very long, for example, 2 miles or more, and is 25 normally composed of a number of individual sections 14 connected end to end. The cable can be straight, slanted or curved. [0014] It will be understood that for simplicity of the drawing a 2D system is shown, whereas in real operations, multiple streamers are generally towed behind an acquisition boat to form a 3D system. Typical 3D vessels tow a suite of such cables, which are separated from each 30 other in the cross-line direction, to allow areal sampling of the wavefield. [0015] The forward section of cable 12 is typically connected to vessel 10 by a lead-in section 16. As is known in the art, each section 12 typically contains a number of hydrophones (not WO 2013/116195 PCT/US2013/023583 shown) that are positioned in each of sections 14 so that they are spaced apart and interspersed. Acoustic wave energy is provided in the vicinity of cable 12 by an air gun 18 or other suitable means. This wave energy travels downwardly through the earth with a portion of it being reflected upwardly at levels where there is a contrast in the acoustic impedance 5 between layers of the strata, for example, at point 20. The hydrophones sense the acoustic pressure waves produced in the water by the upwardly traveling seismic reflections. The propagating wave front also causes movement of the water particles as it passes; the acceleration of the water particles is sensed by the accelerometers. [0016] Referring now to Figure 2, geometry of the ghosting effect is illustrated. In Figure 2, a 10 receiver (e.g. a hydrophone) 32 is located a distance z, below the sea surface 30 and receives a plane wave 34 arriving from below. Wave 34 is travelling at an angle 0 with respect to horizontal. Arrow 36 indicates the primary signal arriving directly at receiver 32. Lines 40, 42, 44, indicate the wave being reflected at the sea surface and traveling down to the receiver, resulting in a ghost signal. An imaginary reflection of the downgoing leg 44 of the ghost 15 signal gives a reflected receiver position 50. From the figure, and using geometry, it can be derived that the signal path for the ghost will be longer than the primary signal path by an amount 6L given by 2z, cos(6) (1) where z, is the receiver depth and 0 is the angle of propagation. 20 [0017] The reflection coefficient at the sea surface is close to -1, so the net effect of the interference of the primary and ghost signals can be represented by the following filter (in the frequency domain) G(co) =1I- exp(ioL /vv) =1I- exp(2icoz, cos(6) /vv) (2) where v, is the propagation velocity in water and a) is the angular frequency. As can be seen, 25 the ghosting filter depends on the angle of propagation. [0018] In preferred embodiments, a "re-datuming filter" is also applied, which relates the data which would be measured by ghost-free receiver at sea level to the ghost-free data measured at depth. The combined filter is given by H (c) = 2i sin(icoz, cos() /v) (3) 30 where i is the unit imaginary number.
    WO 2013/116195 PCT/US2013/023583 [0019] According to preferred embodiments of the present method data recorded by a single receiver, i.e. a "trace," is replaced with a suite of artificial "replacement traces" at nearby locations, which represent the recorded wavefield. [0020] One embodiment of the process is illustrated schematically in Figures 3 and 4. 5 Starting with the initial data in panel 60, a trace selected as the input trace. By way of example only, in Figure 3 and 4, the trace at position 2 is selected as the input trace. The input trace is used to calculate a suite of n angle traces 62, each having with a predetermined time dip (px,py)= (dT/dx, dT/dy). Each time dip corresponds to a distinct propagation angle 0. In 3D, dips in both horizontal directions are sampled. 10 sin(O)= -v p + p (4) [0021] Because the propagation angle of each dip component is known, its ghosting/re datuming filter H can be obtained from Equation (3), and a suitably stabilized inverse can be applied to undo the distortion. Thus, each angle trace 62 can be deghosted with an appropriate, specific deghosting filter, yielding a set of n deghosted angle traces 63. 15 [0022] Each dip component (deghosted angle trace) is then "copied" to selected nearby output locations by applying its time dip, as illustrated by the "Output R positions" 64. In 3D, the output locations are preferably areally distributed around the receiver. The output locations may, but do not necessarily, correspond to the input locations. If the number of selected output locations is also equal to m, this step yields a set of m x n corrected angle traces 64. 20 [0023] The corrected angle traces 64 at each output location are added to give m replacement traces 66. Finally, the m replacement traces 66 from all of the input traces are summed, to give a output trace 68 for each output location. Regardless of the propagation angle of the wave, the output traces 68 will have less signal distortion attributable to ghosts than is present in the original input traces. 25 [0024] In summary, the steps of the preferred method are as follows: a) for each input trace, copy it over to a number of angle traces (the number is preferably determined by sampling theory); b) correct each angle trace for ghosting by using the inverse of the filter H, using Eq. 4 to determine angle 0 and then Eq. 3 to determine the filter 30 H; c) "copy" each deghosted angle trace over to nearby output trace locations, applying its time dip; WO 2013/116195 PCT/US2013/023583 d) at each output location, accumulate contributions from all of the angle/input-trace combinations. [0025] The method described herein, which can be referred to as "beam deghosting," effectively removes the angle-dependent signal distortion regardless of the propagation angle 5 of the wave. In particular, the present beam deghosting technique provides at least three unique advantages: 1) It does not require that the receivers be regularly distributed in space. 2) It can be applied to acquisition systems in which the depth of the receivers varies or not along the cable. The process is input-oriented, i.e. the depth used in Eqs. 1 and 2 is the 10 depth of the input receiver whose contribution is being calculated, thus the procedure naturally accommodates variable receiver depths. 3) It is fully 3-dimensional, i.e. it accounts for waves whose propagation direction has a horizontal component along the direction orthogonal to the cable, as well as a horizontal component in the direction of the cable. 15 [0026] The combination of re-datuming and de-ghosting ensures that the method can be used in systems in which the depth of the receivers varies along the cable because the replacement traces are referred to a common effective receiver depth (zero) before they are summed. [0027] While the foregoing text describes output traces that are on an (x,y) grid, it will be understood that the input traces can however be at arbitrary locations. 20 [0028] In addition, the present methods can be used as a deghosting method for source ghosts. As above, the methods may be applied in 2D or 3D, depending on the sampling of the source locations. [0029] The invention has been described herein with respect to preferred embodiments. It will be understood that modifications and variations can be made without departing from the 25 scope of the invention, which is set out in the claims that follow. 16
    权利要求:
    Claims (6)
    [1] 1. A method for generating a seismic image, comprising the steps of: a) providing at least one input trace corresponding to seismic signals received at a receiver at a known receiver depth; b) using the input trace to generate a series of angle traces, each having a known time dip (px, py)= (dT/dx, dT/dy); 5 c) for each angle trace, using sin(6) = 1vp, + p' where v, is the propagation velocity in water to determine an angle 0 and H(o)= 2isin(icoz, cos(6)/vv) where z, is the receiver depth and ow is the angular frequency to determine a filter H; d) for each angle trace, correcting the data for ghosting by using the inverse of the 10 filter H so as to generate a plurality of deghosted angle traces; e) copying each deghosted angle trace over to a plurality of output locations, applying in each case the time dip for that angle trace, so as to generate a plurality of corrected angle traces; f) at each output location, summing a plurality of the corrected angle traces at 15 that output location so as to generate a replacement trace; and g) using the replacement trace to generate a seismic image.
    [2] 2. The method according to claim 1, further including repeating steps a) - g) for a plurality of input traces, wherein step g) comprises 20 g1) at each output location, summing the replacement traces at that output location so as to generate an output trace; and g2) using the output traces to generate a seismic image.
    [3] 3. The method according to claim 2 wherein the output locations correspond to the input locations. 7 WO 2013/116195 PCT/US2013/023583
    [4] 4. The method according to claim 2, further including repeating steps a) - i) after a period of time so as to generate a second seismic image and comparing the seismic images so as to generate a time-lapse image.
    [5] 5. The method according to claim 2 wherein the number of input traces is greater than
    [6] 100. 6. The method according to claim 1 wherein step b) includes generating at least 30 angle traces from each input trace. 7. The method according to claim 1 wherein the seismic signals are received at a receiver selected from the group consisting of hydrophones, geophones, fiber optic sensors, ocean bottom sensors, streamers, and combinations thereof. 8. The method according to claim 1 wherein the seismic signals are received at a towed marine receiver.
    类似技术:
    公开号 | 公开日 | 专利标题
    AU2013215333B2|2015-01-22|Angle-dependent deghosting of seismic data
    AU2005282945B2|2010-02-18|System for the attenuation of water bottom multiples in seismic data recorded by pressure sensors and particle motion sensors
    US7675812B2|2010-03-09|Method for attenuation of multiple reflections in seismic data
    Wang et al.2012|Premigration deghosting for marine towed streamer data using a bootstrap approach
    US10088588B2|2018-10-02|Device and method for stable least-squares reverse time migration
    GB2503325A|2013-12-25|Receiver-based ghost filter generation
    CA2780081A1|2011-12-15|Method and apparatus for deghosting seismic data
    EP2187240B1|2011-10-26|Method for Optimum Combination of Pressure and Particle Motion Sensors for a 3-D Spread of Dual-Sensor Marine Seismic Streamers
    WO2016076917A1|2016-05-19|Creating a high-resolution earth model using seismic tomography and impedance inversion
    EP2726907A1|2014-05-07|Attenuating sea-surface ghost wave effects in seismic data
    GB2425839A|2006-11-08|Deghosting marine seismic data from sensors at two water depths
    WO2014195508A2|2014-12-11|Systems and methods for de-noising seismic data
    AU2014343360B2|2016-05-26|Hybrid deblending method and apparatus
    EP2491430A2|2012-08-29|Methods to process seismic data contaminated by coherent energy radiated from more than one source
    EP2598917A2|2013-06-05|Methods and systems to eliminate undesirable variations in time-lapse seismic surveys
    CA2815993A1|2013-11-18|Method and apparatus for pre-stack deghosting of seismic data
    WO2014195467A2|2014-12-11|Method and system for simultaneous acquisition of pressure and pressure derivative data with ghost diversity
    CA2827070A1|2014-03-19|Interference noise attenuation method and apparatus
    Cambois et al.2009|Dual-sensor streamer data: Calibration, acquisition QC and attenuation of seismic interferences and other noises
    Amundsen et al.2017|Broadband seismic over/under sources and their designature-deghosting
    AU2014204760A1|2015-07-16|Premigration deghosting for marine streamer data using a bootstrap approach in Tau-P domain
    WO2016128786A1|2016-08-18|Seismic data processing including variable water velocity estimation and compensation therefor
    AU2014201114B2|2014-07-31|Method and apparatus for deghosting seismic data
    AU2014200095A1|2014-07-31|Ghost compensation in beam migration
    Whitmore2016|Extending the 3D primary image with multiples and mirrors
    同族专利:
    公开号 | 公开日
    US9256925B2|2016-02-09|
    GB2512535A|2014-10-01|
    CA2862592A1|2013-08-08|
    GB2512535B|2017-08-02|
    GB201412106D0|2014-08-20|
    WO2013116195A1|2013-08-08|
    AU2013215333B2|2015-01-22|
    US20150003689A1|2015-01-01|
    引用文献:
    公开号 | 申请日 | 公开日 | 申请人 | 专利标题
    US5774417A|1996-10-25|1998-06-30|Atlantic Richfield Company|Amplitude and phase compensation in dual-sensor ocean bottom cable seismic data processing|
    US6539308B2|1999-06-25|2003-03-25|Input/Output Inc.|Dual sensor signal processing method for on-bottom cable seismic|
    GB2397907B|2003-01-30|2006-05-24|Westerngeco Seismic Holdings|Directional de-signature for seismic signals|
    US7450468B2|2004-10-19|2008-11-11|Landmark Graphics Corporation|Methods, systems, and computer readable media for interactively determining parameters for ocean-bottom hydrophone/geophone receiver combination deghosting|
    US7584056B2|2005-02-22|2009-09-01|Paradigm Geophysical Ltd.|Multiple suppression in angle domain time and depth migration|
    US8477561B2|2005-04-26|2013-07-02|Westerngeco L.L.C.|Seismic streamer system and method|
    US8120991B2|2006-11-03|2012-02-21|Paradigm Geophysical S.A.R.L.|System and method for full azimuth angle domain imaging in reduced dimensional coordinate systems|
    FR2955396B1|2010-01-15|2013-03-01|Cggveritas Services Sa|DEVICE FOR PROCESSING SEISMIC MARINE DATA|
    FR2961316A1|2010-06-10|2011-12-16|Cggveritas Services Sa|PROCESS FOR PROCESSING SEISMIC MARINE DATA|
    US9103942B2|2011-10-28|2015-08-11|Westerngeco L.L.C.|Methods and systems for survey designs|
    US9234977B2|2012-01-12|2016-01-12|Westergeco L.L.C.|Processing collected survey data|US9274239B2|2012-01-13|2016-03-01|Westerngeco L.L.C.|Wavefield deghosting|
    US10247822B2|2013-03-14|2019-04-02|Navico Holding As|Sonar transducer assembly|
    US9678235B2|2013-07-01|2017-06-13|Pgs Geophysical As|Variable depth multicomponent sensor streamer|
    US10597130B2|2015-01-15|2020-03-24|Navico Holding As|Trolling motor with a transducer array|
    US11209543B2|2015-01-15|2021-12-28|Navico Holding As|Sonar transducer having electromagnetic shielding|
    US10018719B2|2015-03-05|2018-07-10|Navico Holding As|Systems and associated methods for producing a 3D sonar image|
    WO2016154104A1|2015-03-20|2016-09-29|Schlumberger Technology Corporation|Single streamer deghosting with extended model space|
    US9977143B2|2015-09-04|2018-05-22|Schlumberger Technology Corporation|Method for attenuating multiples|
    US10719077B2|2016-10-13|2020-07-21|Navico Holding As|Castable sonar devices and operations in a marine environment|
    法律状态:
    2015-05-21| FGA| Letters patent sealed or granted (standard patent)|
    2020-08-27| MK14| Patent ceased section 143(a) (annual fees not paid) or expired|
    优先权:
    申请号 | 申请日 | 专利标题
    US201261593090P| true| 2012-01-31|2012-01-31||
    US61/593,090||2012-01-31||
    PCT/US2013/023583|WO2013116195A1|2012-01-31|2013-01-29|Angle-dependent deghosting of seismic data|
    [返回顶部]